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Abstract 

We consider gauge invariant quark two-point Green's functions in which the gluonic phase 
factor follows a skew-polygonal line. Using a particular representation for the quark 
propagator in the presence of an external gluon field, functional relations between Green's 
functions with different numbers of segments of the polygonal lines are established. An 
integral equation is obtained for the Green's function having a phase factor along a single 
straight line. The related kernels involve Wilson loops with skew-polygonal contours and 
with functional derivatives along the sides of the contours. 
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The objective of this work is to investigate the possibilities of deriving integral or 
integrodifferential equations for the two-point gauge invariant quark Green's function 
(GIQGF). Gauge invariant objects seem to have a firmer basis to describe non-perturba- 
tive properties of QCD, as compared to gauge variant ones [HE]- A typical example of 
those is the Wilson loop [3J, which allows the formulation of the confinement property of 
quarks with a relatively simple criterion [HE]. 

Gauge invariant Green's functions involve path-ordered gluon field phase factors. Here, 
we concentrate on two-point quark Green's functions, in which the gluonic phase factor 
follows in general a skew-polygonal line. 

The starting point of our investigation is a particular representation for the quark 
propagator in the presence of an external gluon field, where it is expressed as a series 
of terms involving phase factors along successive straight lines. Then the corresponding 
quantized Green's function becomes expressed in terms of Wilson loops having skew- 
polygonal contours [5J. 

We begin by introducing definitions and conventions. We consider a path-ordered 
gluon field phase factor along a line C yx joining a point x to a point y, with an orientation 
defined from x to y: 



-ig [ V dzTApiz) 

J X 



U(C yx ;y,x) = U(y,x)=Pe J* " w . (1) 

Parametrizing the line C with a parameter X, C = {x(A)}, < A < 1, such that x(0) = x 
and x(l) = y, a variation of C induces the following variation of U [U(x(X),x(X')) = 
U(\,\%A(x(\)) = A(\)] P: 

5U(1,0) = -ig5x a (l)A a (l)U(l,0) + igU(l,0)A a (0)5x a (0) 

+ig f 1 dXU(l, A)x ,/3 (A)F« a (A)5x a (A)f/(A,0), (2) 
J o 

where x' — |f and F is the field strength, F^ = d^A u — d u A^ + iglA^, A v \. 

For paths defined along rigid lines, the variations inside the integral are related, with 
appropriate weight factors, to those of the end points. Considering now a rigid straight 
line between x and y, an ordinary derivation at the end points yields: 

dU Q V ya X) = -tgA a (y)U(y, x) + ig(y - xf £ dX A 17(1, X)Fp a (X)U(X, 0), (3) 
- +igU(y, x)A a {x) + %g{y - xf f~ d\ (1 - A) U(l, X)F Pa (X)U(X, 0). (4) 



dx a 

We adopt the following conventions to represent the contributions of the internal parts 
of the integrals: 

= ig{y _ x f £ dX \ U(l, X)F pa (X)U(X, 0), (5) 
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= ig(y - xf jf 1 dA (1 - A) 17(1, A)F /Jo (A)Z7(A, 0). (6) 

The Wilson loop, denoted $(C), is denned as the trace, in color space, of the phase 
factor (JTJ along a closed contour C: 

$(C) = ^trPe ^ M ; . (7) 
Its vacuum expectation value is denoted W{C): 

W(C) = <$(C)>, (8) 

the averaging being defined in the path integral formalism. The properties of the Wilson 
loop were studied in a long series of papers [Zl El ISl IIOl II^l [T31 1131 IISl IISl IIZl IIHl [13 - 
We adopt the following functional representation for W(C): 

W(C) = e F ( C \ (9) 

In perturbation theory, F(C) is given by the sum of all connected diagrams, the connection 
being defined with respect to the contour C [16]. For large contours and large N c , F(C) 
is proportional to the area of the minimal surface with contour C [H [9] . 

If the contour C is a skew-polygon C n with n sides and n successive marked points 
x±, X2, . . ., x n at the cusps, then we write: 

W( Xn , X n -!, ...,X!)=W n = • • • = e F -. (10) 

The two-point GIQGF, with a phase factor along a line C, is defined as 

S a/3 (x,x';C x , x ) = -— ($ p (x')U(C x > x ;x\x) i) a (x)) . (11) 

(The color indices are implicitly summed.) For skew-polygonal lines with n sides and 
n — 1 junction points yi, 7/2, • • -,2/n-i between the segments, we define: 

S( n )(x,x';y n -!,...,y!) = -— (^(x')U(x', y n -i)U(y n -i, y n -2) ■ ■ ■ U(y u x)tp(x)). (12) 

For one straight line, one has: 

S (1) (x,x') = S(x,x') = ~ ®(x')U(x',x)4>(x)). (13) 

( The index 1 will generally be omitted from that function.) 

We shall adopt a two-step quantization method. One first integrates with respect to 
the quark fields. This produces in various terms the quark propagator in the presence of 
the gluon field. Then one integrates with respect to the gluon field through Wilson loops. 
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To make Wilson loops appear, one needs an appropriate representation for the quark 
propagator in extenal field. We use the following representation which involves phase 
factors along straight lines together with the full quark Green's function = S 0, [21] : 

where S(x, x'\ A) is the quark propagator in the presence of the external gluon field A. 
(Integrations on intermediate points are implicit.) This equation yields an expansion of 
S(x, x'; A) in terms of the GIQGF S and explicit phase factors along straight lines. The 
above representation is a generalization of that used for heavy quark propagators starting 
from the static case [22] . 

Systematic use of Eq. ffl~4")) leads to functional relations between various GIQGFs. 
Consider for this the Green's function 5( n ) [Eq. (fl2l) ]. Integrate with respect to the quark 
fields; one obtains: 

S(n) 0, x'; y n -i, . . . , yi) = — (U (x', y n - 1 )U(y n - 1 , y n ^ 2 ) ■■■V {yi, x)S{x, x'; A)). (15) 
Use of the expansion (TBI) for S(A) gives: 

S {n) (x, x'- y n _ u . . . , yi ) = S(x, x') e F «+i( x '' Vn-u ■■■^ x ) 

+ { ij a+ + S{x,y n )—z)l S {n+1 )(y n ,x;y n -i,...,y 1 ,x). (16) 



Notice the appearance of the Wilson loop average along the skew-polygonal contour with 
(n + 1) sides. Repeating the same operation for S^n+i) and so forth, one can in principle 
express any St n ) in terms of the lowest-order Green's function S and Wilson loop averages 
and their derivatives along the corresponding contours. 

Next, we use the equation of motion of the quark fields. This yields for S( n y. 

(ij.d {x) - m)S {n) (x, x'; y n _!, . . . , Vl ) = i5\x - x')e F ^ x > f»-i> 

Sx^ ' 1 ' 

We observe that the right-hand side contains, as an unknown, the functional derivative 
along the rigid segment xy\ of Si n ). The task of evaluating that quantity is, however, 
facilitated by the functional relations (fTBl . which relate two successive Green's functions 
with increasing index. They allow the evaluation of the rigid path derivative of a Green's 
function in terms of a similar derivative of a Wilson loop average and the derivative of a 
Green's function with a higher index. Systematic repetition of this procedure allows us 
to express the derivative of a Green's function in terms of a series of Green's functions 
whose coefficients are functional derivatives of Wilson loop averages. One thus obtains 
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chains of coupled integral (or integrodifferential) equations between the various Green's 
functions. At the end, each Green's function S/ n ) can be expressed, at leading order of an 
expansion, by means of the functional relation (Tl6|) . in terms of the lowest-order Green's 
function S. Thus an equation where solely the Green's function S would appear becomes 
reachable. 

In the present work we are mainly interested by the simplest Green's function S; the 
general structure of the derivative 8S/8x^~ is: 

^ '_ = K 1 ^(x',x)S(x,x') + K 2 ^(x / ,x,y 1 )S {2) (y 1 ,x';x) 

oo 

+ K n^-{x\x,y u . . .,y n -i) S( n )(y n - h x';x,yi, . . .,y n -2)- (18) 

n=3 

The kernels K n (n = l,...,oo) are composed of Wilson loop averages along skew- 
polygonal contours with (n + 1) sides, of their derivatives along these sides and of (n — 1) 
quark Green's functions S and their derivatives. The total number of derivatives con- 
tained in K n is n. Each segment of the contour is submitted at most to one derivation. 
The calculation of the expression of the kernel K n requires solely consideration of terms 
of order lower or equal to n. 

The explicit expression of 5S/5x tl ~ up to the third-order of its expansion is: 

SS(x,x') 5F 2 {x',x) , ^F^x'.x.yx) 

= S{x, x ) - -^_- yr S(x, Vl ) 7 S (2) ( yi , x ; x) 

5 3 F 4 (x',x,y 1 ,y 2 
5x"-5y^ 1+ 5y^ 

+ ••• . (19) 

The contents of the kernels K n can also be classified according to well-known structures, 
such as completely connected, crossed, nested, etc. 

The integral form of the equation of motion (1171) for n — 1 is: 

S(x,x') = S (x,x') + J d*x"S (x,x")r 5S ^'\ (20) 

in which one has to inject the expression of SS/Sx~ resulting from Eq. ( Tl8l) . (So is the 
free quark propagator.) 

At short distances, governed by perturbation theory, each derivation introduces a new 
power of the coupling constant and therefore the dominant terms in the expansion are the 
lowest-order ones. At large distances, Wilson loops are saturated by the minimal surfaces 
having as supports the contours [SI [HI US, HI EDI EI]- Here also, the dominant contributions 



+ ^^t^, ai ^]]a 2 2 + S ( x ' Vi) 7 Ql S{Vi, 2/2) 7° 2 S , (3)(j/2 > x, yi) 
S 2 F 4 
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come from the lowest-order derivative terms. Therefore, the expansion in Eq. (1201) can be 
considered in general as perturbative whatever the distances are, provided that for each 
type of region the appropriate expressions are used for the Wilson loops. The first term 
of the expansion, represented by a single derivative, is null for symmetry reasons (the 
derivative SF^x', x)/8x~ being orthogonal to xx'). Hence the non-zero dominant term 
of the expansion is the second-order derivative one. Furthermore, the various Green's 
functions S( n ) are themselves dominated by their lowest-order expression of Eq. ffTTl) . 
involving only S and a Wilson loop. In that approximation, 5S(x, x')/5x~ takes the form 

Thus, the dominant part of the integral equation (120]) to be solved reduces to a closed 
form relative to the full Green's function S. 

The approximate form (j2ip can be used for a first resolution of the integral equation 
(GUI). 
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